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We consider critical models in one dimension. We study the ground state in thermodynamic limit 
[infinite lattice]. Following Bennett, Bernstein, Popescu, and Schumacher, we use the entropy of a 
sub-system as a measure of entanglement. We calculate the entropy of a part of the ground state. 
At zero temperature it describes entanglement of this part with the rest of the ground state. We 
obtain an explicit formula for the entropy of the subsystem at low temperature. At zero temperature 
we reproduce a logarithmic formula of Holzhey, Larsen and Wilczek. Our derivation is based on 
the second law of thermodynamics. The entropy of a subsystem is calculated explicitly for Bose gas 
with delta interaction, the Hubbard model and spin chains with arbitrary value of spin. 
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Entanglement is an im por tant resource f or q uantum 
computation^ SIllilElliiliilsllslillslIll 
l59| . Recently entanglement was studied in 
detail in ground state of different physical models 0, Is^ 

m S Ei E m El . 

Conformal field theory is useful for the description 
of low temperature behavior of gap-less models in one 
space and one time dimensions. Critical models are clas- 
sified by a central charge c of corresponding Virasoro 
algebra. A definition of the central charge is given the 
the Appendix. Conformal field theory is closely related 
to Luttinger liquid 8] approach. We are interested in 
specific entropy s [entropy per unit length] . Let us start 
with specific heat C = Tds/dT. Low temperature be- 
havior was obtained in [(A 

Here c is a central charge and v is Fermi velocity. We 
put both Plank and Boltzmann constants equal to ks = 
h — 1. Later in the paper we shall use: 



C 



ttTc 



as T ^ 0. 



(1) 



We are more interested in specific entropy s. We can 
integrate the equation and fix the integration constant 
from the third law of thermodynamics (s = at T = 0). 
So for specific entropy we have the same low temperature 
behavior: 



TiTc 



T -> 0. 



(2) 



For quantum spin chains this formula agrees with [ill ] . 
To formulate the problem precisely let us consider Bose 



gas with delta interaction. The Hamiltonian of the model 



is: 



H = 



do; 

0. 



(3) 



Here t/j is a canonical Bose field and q > is a coupling 
constant. The model was solved in [13]. Physics of the 
model is described in a book First let us consider 

the model at zero temperature and in the infinite volume. 
The ground state is unique \gs). We consider gas with 
positive density. We are interested in the entropy S{x) 
of the part of the gas present on a space interval (0,x). 
Formally we can define it by means of the density matrix 



P = troo {\gs){gs\) 



(4) 



Here we traced out the 'external' degrees of freedom, they 
describe the gas on the rest of the ground state: on the 
unification of the intervals (—00, 0) and (cc, 00). The den- 
sity matrix p describes gas on the interval (0, x). Now we 
can calculate von Neumann entropy S{x) of the part of 
the gas on the interval (0, x) : 



S{x) — — tropin p 



(5) 



Here we are taking the trace with respect to the degrees 
of freedom representing the part of the gas on the in- 
terval (0,a;). In the major text books it is shown that 
the laws of thermodynamics can be derived from statis- 
tical mechanics, see for example 0, 0] . Second law of 
thermodynamics states that the entropy is extensive pa- 
rameter: the entropy of a subsystem S{x) is proportional 
to the system size x: 



Six) 



at T > 



(6) 
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Thermodynamics is applicable to the subsystem of 
macroscopical size, meaning large x. Here specific en- 
tropy s depends on the temperature. For small temper- 
atures the dependence simpHfies, see (2): 

S{x) — — — X, X > —. (7) 
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Let us try to find out how S{x) depends on x for zero 
temperature. It is some functions of the size: 



S{x) = f{x), at r = 0. 



(8) 



Now let us apply the ideas of conformal field theory, see 
Hi, [13 and also Chapter XVIII of (IJl . We can arrive 
to small temperatures from zero temperature by confor- 
mal mapping exTp[2'!TTz/v]. It maps the whole complex 
plane of z without the origin to a strip of the width 1/T. 
This replaces zero temperature by positive temperature 
T. The conformal mapping results in a replacement of 
variable x by [v/nT] sinh[7rTa;/f]. So the entropy of the 
subsystem at temperature T is given by the formula: 



S{^) = / ( ^ sinh 



ttTx 



at T > 



(9) 



In order to match this to formula (7) we have to consider 
asymptotic of large x. The formula simplifies: 



Six) ^ f i exp 



nT{x — Xq) 



Tx 



oo. 



(10) 



Here ttTxq/v = -ln(u/27rr) . 

Formulae (7) and (10) should coincide. Both represent 
the entropy of the subsystem for small positive tempera- 
tures. This provides an equation for /: 



/ exp 



7rT(x — Xq) 



^(^-^o) (11) 



This formula describes asymptotic for large so we 
added — xq to the right hand side. We are considering 
the region x > 1/T and xq ^ ln(l/r), so x >> xo at 
T — > 0. In order to solve the equation for /, let us intro- 
duce a new variable y — exp [ttT(x — Xq)/v] . Then the 
equation (11) reads: 



/(y) = o 



(12) 



So we found the function / in (8). Now we know that 
at zero temperature entropy of the gas containing on the 
interval (0, x) is: 



S (x) = 2 Ina:^, 



as X 



(13) 



Let us remind that for Bose gas the central charge c — 1, 
see 10 . Our result agrees with the third law of thermo- 
dynamics. Specific entropy is a limit of the ratio S{x)/x 
as a; — > oo. The limit is zero. 

Now we can go back to our formula © and substitute 
the expression for /, which we found: 



c ( V 

s{x) = 3 In i:;^^^"^^ 



irTx 



(14) 



This formula describes crossover between zero and large 
temperature. 



The proof presented here is universal. It is also appli- 
cable to quantum spin chains. For example to XXZ spin 
chain. The Hamiltonian of the model is: 



y 

i+i 



2/icr|} (15) 



It describes interaction of spins 1/2. Here a are Pauli 
matrices. The model can be described by a conformal 
field theory with the central charge c = 1 in a disordered 
regime: — 1 < A < 1 and magnetic field h is smaller 
then critical he, see jl^. For magnetic field larger then 
critical the ground state is a product of wave functions 
of individual spins. Entropy of any subsystem is zero. 

The formula (13) was discovered in "^,51. More general 
formula can be found in (36n37| |. The entropy of the sub- 
system S{x) describes the entanglement of the subsystem 
and the rest of the ground state. 

The formula H13(l describes the entropy of large sub- 
system X — > oo at zero temperature, first correction to 
this formula is a constant term: 

S (x) = - In — , as a; ^ oo 

3 Xq 

The constant term defines the scale xq. It does not de- 
pend on x, but it does depend on magnetic field and 
anisotropy, see [l^ . 

So far we discusses the entropy of entanglement for 
critical models. For non-critical [gap-full] models asymp- 
totic of S{x) is a constant, depending on parameters of 
the model. For AKLT spin chain S'(oo) ~ 2, see 01 • For 
XY spin chain S{oo) depends on anisotropy and mag- 
netic field and shows singularities at phase transitions 
see mill. 

We evaluated von Neumann entropy (0). Calculation 
of generalized entropies [Renyi and Tsallis] can be re- 
duced to calculation of a trace of some power a of the 
density matrix: tvxip") it was done in the paper |l6|. 



SPIN CHAINS WITH ARBITRARY SPIN 

In this section we consider higher values of spins. We 
calculate dependence of entropy of a subsystem on the 
value of interacti ng sp in s. Let us consider isotropic XXX 
anti-ferromagnet [23. For spin s = 1/2 we can represent 
the Hamiltonian as: 

Hi=E„^n, (16) 

= SnSn+l = Sf^Sf^j^i + -|- S^S^^i (17) 

The model is solvable by Bethe Ansatz l^^. A general- 
ization of this model to spin s = 1 was found in 



Hi - ^ [Xn - Xl] 



(18) 
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It was solved by Takhatajan and Babujian, see 
Generalization for higher spin s was found in 




(19) 



The function F{X) is a polynomial of a degree 2s. It can 
be written as 



Fix) 
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1 

I n 



X - Vi 

^ ^ k yi — Vj 

' =0 k=l+l j ^ 



(20) 



Here yi = 1(1 + l)/2 - s{s + 1). The model lO also solv- 
able by Bcthe Ansatz. There is no gap in the spectrum 
of this Hamiltonians. The models can be described by a 
conformal field theory with the central charge: 

3s 



s+1' 

see [3^. So the entropy of a block of x spins is: 
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Inx, as X 



00 



(21) 



The coefficient in- front of the In a; increases from 1/3 to 
1 as spin s increases from 1/2 to 00. 

Note. It is interesting to remark that in ferromagnetic 
case [different sign in from of the Hamiltonian] the en- 
tropy also scales logarithmically, but the coefficient in 
front of the \nx is different, see 



see Also Fermi velocities are different for spin Vs 

and charge Vc degrees of freedom. Both spin and charge 
degrees of freedom contribute to the specific entropy s. 
So for small temperature we have: s = ttT/Svs + ttT/Svc- 
The entropy of the subsystem is proportional to the size 
of the system (second law): 



(23) 



Now we have to apply conformal arguments separately 
to spin and charge degrees of freedom. The results add 
up: 
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-\nx 



(24) 



This describes the entropy of electrons on the interval 
(0,a;) in the infinite ground state at zero temperature. 
This is actually an asymptotic for large x. In case of 
the Hubbard model crossover formula [which mediate be- 
tween zero and positive temperature] looks like this: 
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(25) 



i In ( sinh 
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i In ^ Sinh 
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The approach developed here is universal, see ^ . It 
is applicable to other models of strongly correlated elec- 
trons, see For example formulae (|24|l . (125(1 describe 
entropy in t-J model below half-filling as well. 



HUBBARD MODEL 

The Hamiltonian for the Hubbard model H can be 
represented as 

j j 

Here is a canonical Fermi operator on the lattice 

[operator of creation of an electron] and rij^a — c] a-'^j.a 
is an operator on number of electrons in cite j with spin 
a. Cite summation in the Hamiltonian goes through the 
whole infinite lattice. Coupling constant m > and mag- 
netic field h is below critical. At half filled band only 
spin degree of freedom is gap-less, so at zero temperature 
Six) = (1/3) In a;. Let us consider the model below half 
filling [less then one electron per lattice cite] . The model 
was solved in |0] . Detailed description of physics of the 
model can be found in the text-book 0^. Charge and 
spin separate in the model. The model is gap-less. Both 
charge and spin degrees of freedom can be described by 
Virasoro algebra with central charges Cc = 1, = 1. 



SUMMARY 

In the paper we described universal properties of the 
entropy in one dimensional gap-less models. We studied 
the entropy of a subsystem. We considered scaling of 
the entropy in spin chains, strongly correlated electrons 
and interacting Bose gas. We discovered crossover for- 
mula for entropy of the subsystem, see It describes 
entropy of a large subsystem at low temperature. We 
discovered a dependence of entanglement on spin s for 
antiferromagnet, see ((21|l . 



APPENDIX 

Important characteristic of conformal field theory is 
a central charge. It can be defined by considering an 
energy-momentum tensor T^^i^iz). Here z is complex 
space-time variable z — x -\- ivt and z; is a Fermi ve- 
locity. The component of energy-momentum tensor with 
fi ~ v ^ z is denoted by Tz.z = T. Correlation function 
of this operator has a singularity: (T(z)r(O)) = c/z**, the 
coefficient c is the central charge. 
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